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Abstract 



It is shown that the radiation of a charged Higgs boson off a third- 
generation quark (charged-Higgs-strahlung) provides an important channel 
for the discovery of the charged Higgs at hadron colliders. Equivalently, 
in supersymmetric models with explicit lepton-number (i?-parity) violation, 
sleptons may also be produced in association with quarks (slepton-strahlung). 
Higgs- and slepton-strahlung production cross sections are given for both the 
Tevatron and the LHC. The LHC cross sections imply that heavy O(TeV) 
charged Higgs bosons can be produced via charged-Higgs-strahlung and that 
strahlung production of charged sleptons is possible even for small i?-parity 
violating couplings. The possible discovery of sleptons through this channel 
offers a surprising handle on models of neutrino masses. 
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I. INTRODUCTION 



While the LEP experiments reach their final phases, the upgraded Tevatron (s/s = 
2TeV) at FNAL and the CERN Large Hadron Collider (LHC) (y/s = 14TeV) will provide 
the next front in searches for particles associated with physics beyond the Standard Model 
(SM). In particular, since LEP cannot probe masses significantly above the 100 GeV mark, 
the probe of heavier particles will remain the task of these two hadronic machines. Sufficient 
production cross sections at these colliders, however, may be difficult to achieve if the new 
particles carry only weak charges. Effective couplings of these particles to (initial) gluons 
may be induced at the quantum level, but they suffer from loop suppression and, hence, 
are generically small. Sufficiently large and measurable production cross sections require 
either some enhancement of the radiative gluonic couplings or an alternative production 
mechanism. 

In the case of the Higgs bosons present in supersymmetric models, as well as non- 
super symmetric two-Higgs-doublet models, there exist two potentially large parameters that 
can partially compensate for the otherwise small couplings: the intrinsically large t-quark 
Yukawa coupling and the ratio of vacuum expectation values (vev's) of the two neutral Higgs 
bosons tan/? = i^/Vl, which is constrained from above to be < 60 by the perturbativity 
of Yukawa couplings. They may (i) enhance the radiatively induced gluonic coupling of 
the neutral Higgs, in proportion to the t-quark Yukawa coupling, leading to production via 
gg — > H° 0; (ii) sufficiently increase the rate for the decay of the t-quark into a (light) 
charged Higgs boson, t — > H + b [@-^]; or (Hi) enhance the Higgs-strahlung associate produc- 
tion through the 2 — > 3 partonic processes gg, qq — > qqH°, qq'H ± [|6| pT0|] , through the 2-^2 



ones qg — > qH°, q'H ± ||ll,12], and through the 2 — ► 1 process bb — *> H° [^,10|. (For an 



overview, see, for example, Ref. JT3]].) 

From charge conservation, a radiatively induced gluonic coupling cannot lead, at the 
level of an elementary process, to the production of only one charged Higgs. Hence, one 
needs to consider either the single production from t-quark decays or the production in 
association with quarks - which is referred to as Higgs-strahlung. Production in association 
with gauge bosons WH°/ZH°, W ± H T , and in supersymmetric models, in association with 
squarks, is also possible. The former mechanism, however, leads to subleading cross sections, 
which are difficult to observe p^Jl5|] . The latter p6^T8[ is more model-dependent and will 



not be considered here. The production mechanisms with one single Higgs (and quarks) 
in the final state, or "single production", are kinematically in advantage with respect to 
pair-production mechanisms. The latter includes (a) the Drell-Yan process qq — > H + H , 
which is suppressed by weak couplings, and, at the LHC, by the low quark luminosities 
(relative to the gluon one); (b) effective gluonic couplings, which are now allowed by charge 
invariance, gg — > H + H ||19| . Pair production, therefore, does not allow the discovery of 
the charged Higgs at the Tevatron. At the LHC, it provides a limited discovery reach of the 
charged Higgs in a non-supersymmetric two-Higgs-doublet model fll9fl , but it may become 
more competitive in supersymmetry where additional contributions to the effective gluonic 



couplings arise 20,21 



Our focus here is on the single production of a charged Higgs boson. Single production 
from t-quark decays plays the most important role when it is kinematically allowed, that 
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is for m H ± < m t — m&, and was studied extensively by various authors [§]]. The charged- 
Higgs-strahlung in the 2 — 3 channels gg, qq — *> tH~b encompasses the resonant production 
of a pair ti followed by the decay i — > H~b in the same kinematical region m#± < m t — 
in addition, it provides the intrinsically off-shell associate production of H ± beyond this 
kinematical limit, throughout all possible ranges of m#±. Strahlung from the 2 — > 2 channel 
gb — > tH~ is also possible, since the 6-quark is obtained from the proton via a gluon. 
Therefore, both types of partonic processes, the 2^3 and the 2 — ► 2, give rise to inclusive 
processes that are formally of the same order in an a s expansion. Note that away from the 
resonance region, H~ decays dominantly into ib, and the i-quark into W + b. Thus, the 2^3 
processes gg, qq — > tH~b give rise to four 6-quarks in the final state, whereas only three 6's 
are produced in the 2 — > 2 process gb — > tH~ . Both elementary processes contribute to 
the inclusive production of a single charged Higgs at hadron colliders, when at most three 
6-quarks are tagged and used to identify the final particle configuration. If four 6's can 
be detected, it is possible to measure each cross section separately. Otherwise, the two 
elementary processes have to be properly combined into an inclusive cross section, avoiding 
double counting of the contributions coming from gb — > tH~ and from gg — > tH~b, when 
one of the two gluons produces a bb pair collinear to the initial proton (or antiproton). Since 
identification and detection issues will not be discussed here, both the individual- channel 
and inclusive cross sections will be given below. 

The charged-Higgs production cross sections are calculated and illustrated for the up- 
graded Tevatron and the LHC in Section |J, where all relevant issues are discussed in detail. 
Attention is given to the prospects of discovery at the Tevatron, which may constrain m#± 
beyond the kinematical limit < m t — m^. On the LHC front, it is found that a charged 
Higgs as heavy as 0(TeV) may be produced at the LHC via the strahlung processes. The 
results shown are valid for the charged Higgs of supersymmetric and non-supersymmetric 
two-Higgs-doublet models. The charged-Higgs decay modes, however, may differ in the two 
classes of models if m H ± is sufficiently large 0.^] . 

For reference, the indirect lower limit on m H ± coming from the measurement of the 



inclusive decay b — > 37 amounts (at present) only to ~ 165 GeV [f22] in non-supersymmetric 
two-Higgs-doublet models. No substantial limit exists for charged Higgs in supersymmetric 
models, when supersymmetric partners can be exchanged in the loop mediating the b — > 57 
decay. Direct lower bounds on m#± are given by collider searches at LEP II and at the 
Tevatron. The LEP II bound, m H ± > 54 GeV ^ at y/s = 130 GeV (m H ± > 68 GeV at 
higher energy runs |13[), applies to the case of a charged Higgs boson present in two-Higgs- 
doublet models. The Tevatron searches give combined (and currently modest) bounds in 
the m#±-tan/3 plane for supersymmetric and non-supersymmetric charged Higgs bosons. 
Both searches have been critically discussed in Ref . || . Conservatively, all results presented 
in this paper are shown for m#± ^ 45 GeV, the model-independent limit extracted from the 
measurement of the Z width. 

In supersymmetric models in which i?-parity and lepton number are not conserved, the 
hypercharge Y = — 1 Higgs and slepton fields are not distinguished by any quantum num- 
bers and could interact in a similar way. Thus, in these models, sleptons can be produced 
via slepton-strahlung just as the (charged) Higgs. The relevant Yukawa couplings, i.e. the 
slepton-fermion-fermion couplings, are subject to various low-energy constraints but are 
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otherwise arbitrary, since they do not relate to fermion masses (once the two Higgs doublets 
are identified as those whose neutral components are aligned along the two large vev's). 
They are, however, related to radiative neutrino masses, as explained below, and there- 
fore measurements of lepton-number-violating operators from slepton production provide a 
unique and important window on this class of models for neutrino masses. The associate 
production cross section for the charged sleptons is presented for both the Tevatron and 
the LHC experiments in Section |TJ. It is shown that charged sleptons can be produced in 
abundance for i?-parity-violating couplings as small as 0.01. Relations and possible lessons 
to models of neutrino masses are also demonstrated in Section |TJ. We also comment on the 
case of the neutral sleptons, the sneutrinos u, which is complicated by the presence of the 
gluon fusion channel gg — > v. 

Results and discussions of the potential impact of these strahlung channels on future 
searches are summarized in Section |V|, where we also comment on the possibility that both 
charged-Higgs and slepton-strahlung channels are present. 

All calculations are done at the leading order in QCD. Higher-order corrections may be 



important, as was shown in the case of associate production of the neutral Higgs ||24|| . Their 
inclusion is called upon, but this is left for future study. 



II. CHARGED-HIGGS-STRAHLUNG 



The charged Higgs boson interacts with quarks according to the Lagrangian 

c = { tan ^) Vji dRi + \M^r cot ^) Vj *' llRi dhj 1 H+ + h ' c ' ' ^ 

where the standard notations for the SU(2) coupling g, the up- and down-quarks Ui and 
di of a generation i, and for the Cabibbo-Kobayashi-Maskawa (CKM) matrix V are 
used. Hereafter, all intergenerational mixing terms are neglected, as well as all Yukawa 
couplings other than those for the t- and 6-quarks, which have respectively strength 
ht ^ gmt/(V2M w sin (3) ~ l/sin/3 and h b ~ gm b /(y/2M w cos/3) ~ 0.017 tan/3. (The 
current mass mj = 3 GeV is used hereafter. Note that using the pole mass instead will 
increase hi,, an affect that is important in the large tan/3 regime.) Since all calculations are 
done to leading order, model-dependent radiative corrections to these relations [|],^| are 
also omitted. They could, however, be large and play an important role (at the order in 
perturbation theory at which they must be included) by smearing the tan j3 dependence of 
the if ^production cross sections. 



A. The 2 — ► 3 process: cross-section calculation 



Quarks produced in quark or gluon collisions can radiate a Higgs boson leading to the 
associate production gg, qq — >■ tH~b. Diagrams describing the quark-initiated parton pro- 
cesses are shown in Fig. [I], those describing the gluon-initiated processes in Fig. 0. As 
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FIG. 1. Diagrams contributing to pp(pp) — ► tH (b)X through an elementary quark-initiated 
2^3 process. 





FIG. 2. Diagrams contributing to pp (pp) — ► ti? (£>)-X" through an elementary gluon-initiated 
2^3 process. 

explained in the introduction, if the center-of-mass energy of the relevant collider is suffi- 
ciently large to allow the tagging of two 6-quarks in addition to the if ± -decay products, the 
measurement of the 2^3 production cross section may be possible. 

The production cross section for a pp (pp) collider is obtained as usual by convoluting the 
hard-scattering cross sections of the quark- and gluon-initiated processes with the quark- 
and gluon-distribution functions inside p and/or p: 



a = 7T — I rfPS (9l +92;Pl,P2,P3) X 

2« JT min T Jt X J 

\J2 [q( x iVf)q( T / x iVf) +1^ Q]\M\ 2 qq + g{x,n f )g{T/x,n f )\M\ 



99 



(2) 



where qi (p^ are the four-momenta of the initial partons (final-state particles), the func- 
tions q(x,fif), q(x,fif) and g(x,fif) designate respectively the quark-, antiquark- and gluon- 
density functions with momentum fraction x at the factorization scale fif, and the index 
q in the sum, runs over the five flavors u, d, c, s, b. Finally, s is the hadron center-of-mass 
energy squared, whereas the parton center-of-mass energy squared is indicated, as usual, by 
s = XiX 2 s = ts, with r min = (m H ±+m t +m b ) 2 / s. The integration over dPS(qi+q 2 ] Pi, P2, P3) , 
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an element of phase space of the 3-body final state, can be reduced to four nested integrals 
with bounds explicitly given in appendix [A]. 

When rewriting the third-generation vertex t-b-H ± in ([[]) as i(g/ \/2) V tb (v + 075) +h.c, 
with vector and axial couplings v and a given by: 

1 / m b m t \ 1 / m b m t _\ , nS 
v — — — — tanp + — — - cotp ; a — — — — tanp — — — cotp , (3) 

2 \M W M w H J 2 \M W M w H ) K ' 

the square amplitudes and l-M]^ can be decomposed as: 

\M\%= ^M 2 W ) (^a s ) 2 \V tb \ 2 (v 2 V™ + a 2 A™) 

\M\1=(a^MI\ (4na s ) 2 \V tb \ 2 (v 2 V™ + a 2 A^) . (4) 

In this notation color factors are included in the reduced squared amplitudes V q ^, A qq and 
V", A", while the strong coupling constant is explicitly factored out. The expressions for 
the reduced amplitudes are too lengthy to be given hereQ. 

In the kinematical region m t > m H ± + m b the above cross section could be well approx- 
imated by the much simpler resonant production cross section, as given by the on-shell tt 
production cross section times the branching fraction for the decay t — > H~b. 

B. The 2^3 process: numerical results 

The production cross section a(pp — > tH~bX) for the Tevatron is shown by the solid lines 
in Fig. |asa function of the charged-Higgs mass for three different values of tan f3, tan (3 = 2, 
10, and 50, and in Fig. |asa function of tan/3 for different values of m H + = 100, 140, and 
200 GeV. Notice that these figures show the cross section for single production of H~ only 
and that identical results are obtained for the process a{pp — > tH + bX). All calculations 
are done at the leading order in QCD. The leading-order parton distribution functions 
CTEQ4L |26[ are used in all calculations, and the renormalization (hr) and factorization 
(fj,f) scales are always fixed to the threshold value mt + mn±. The variation of these scales 
results in general in changes in the cross section presented here: a variation in the interval 
between {m t + m H ±)/2 and 2{m t + m H ±) can produce deviations up to ±30% with respect 
to the values shown in the figures. Higher-order corrections, therefore, may be important, as 
was shown in the case of associate production of the neutral Higgs. Their inclusion is called 
upon, but it is left for future work. As a cross-check of our cross-section calculation and 
phase-space integration procedure, we have reproduced, by taking the appropriate limit, the 
well-known gg, qq tih cross section [|] for the neutral Higgs h production in association 
with t-quarks. 



The Fortran code for these amplitudes is available upon request. 
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FIG. 3. The leading-order production cross section a(pp — * t(b)H~X) as a function of the 
charged Higgs mass, for y/s = 2TeV, is shown for three different values of tan/3 = 2, 10, 50. The 
solid lines indicate the cross sections obtained from 2 — ► 3 elementary processes gg — * tbH~ and 
qq — > tbH~~; the dashed lines show the cross sections obtained from the 2 — ► 2 process gb — * tH~ . 
Renormalization and factorization scales are fixed as hr = fj,f = mt + rnjj-. 



t — i — i — rq 




j i i i_d 



10 20 30 40 50 

tan /? 

FIG. 4. The production cross section a(pp — > t(b)H~X) is shown as a function of tan/3 for 
mjj± = 100, 140, and 200 GeV (from top to bottom). As in Fig. p| solid lines correspond to the 
2 — > 3 processes, dashed lines to the 2 — > 2 process. 
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FIG. 5. The leading-order production cross section a(pp — > t(b)H~ X), for y/s = 14TeV is 
shown as a function of the charged Higgs mass for three different values of tan/3 = 2, 10, 50. 
As in Figs. ^ and the solid lines indicate the cross sections obtained from 2 — ► 3 elementary 
processes, the dashed lines, the cross sections obtained from the 2^2 process. Renormalization 
and factorization scales are fixed as fj 1 R=fj 1 f = mt +mjj-. 

As expected, the production cross section is enhanced in the resonance region m#± < 
m t — m b , when H ± is obtained as a decay product of one of two t-quarks produced on 
shell (see first diagram in Fig. [3] and the first and third diagrams in Fig. 0). The resonant 
i-quark propagator is regularized by the width of the t-quark, calculated from the SM decay 
t -> bW + : 



F(t -> bW~ 



64tt 



I V* 1 2 
I v tb\ 



m t M? 



w+ 



M 2 W + 



(m 2 + mf) + (m 2 t -m 2 b y 



A 1/2 1 



«%+ m% 



mi 



(5) 



where A is the Kallen function X(x,y,z) = ((x 2 
t -> H + b: 



v2\2 



4:yz); and from the decay 



r(t -> 



32tt 



m t J 



m 



H+ 



ml 



+ aT 



1 



A 1/2 1 



m 



m 



2 ' 



-4 , (6) 



mi 



for each value of m#± and tan/5. In this region the cross section is not distinguishable 
from the convolution of a(pp — > tiX) with the branching fraction Br(t — ► H + b). The tan/3 
dependence, shown explicitly in Fig. f|, has the same typical pattern as the branching ratio 
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Br(t — > H + b), i.e. large enhancements for very small and very large values of tan/3 and a 



minimum for tan (3 ~ Jmt/mi,. Away from the resonance region, the cross section diminishes 
rather rapidly and becomes negligible (at the Tevatron energy) for m#± ^ 250 GeV. 

Results qualitatively similar to those found for the Tevatron center-of-mass energy are 
obtained in the case of LHC searches. They are shown in Fig. [5[ Assuming integrated 
luminosity of 100 fb _1 , a significant production cross section is obtained for m#± ^ 1 TeV, 
even for tan (3 = 10, which is near the minimum of the production cross section. 

Our results generalize those of Ref. 0, where contributions from the gluon- initiated 
diagrams were calculated for the LHC but only for m t = 150 GeV and for a fixed tan (3 = 1. 
Both of these values are now experimentally ruled out. We also note a disagreement between 
our calculation and that of Ref. 0, which could be, in part, due to our usage of a current set 
of structure functions. A similar calculation for the LHC case was also presented in Ref. [§ , 
for different values of m t . No immediate comparison with our results is, however, possible. 



The mechanism of associate production was also emphasized in Ref. [27], and a preliminary 



study was presented for the Tevatron center-of-mass energy [28 . 



C. The 2^2 process 




FIG. 6. Diagrams contributing to pp (pp) — * tH X through an elementary 2 — ► 2 process. 



As mentioned in the introduction, for an inclusive measurement of single production of 
H~, or if only one 6-quark can be tagged besides the decay products of H~, the elementary 
process gb — > tH~ has to be considered as well. The corresponding diagrams are shown in 
Fig. H and the hard-scattering cross section reads: 



a(gb -> H~t) = ( 4%M^ ] (4na s ) \V t 



V2 



1927TS (l-x 2 b f 1 



L(l-x 2 



2b 



+ 



C 4 



2L (l +xl —2x\x\ t —2x\ t {\- 



~ x \t, 



S i X ht 



b(3-7x 2 ht +xi(3+xl t )+2x 2 b (l-x 2 ht ))]} , (7) 
[(l-(x t + x,) 2 )(l-(x i -x,) 2 )] 1/2 , 



x 2 ; b 



where C± = v 2 ± a 2 ; Xi = m, 
and L = In [(1 - x\ t + 6)/(l - x\ t - b)\. 

The cross section originated from the 2^2 process only is shown in dashed lines in 
Figs. ^| and |] for the upgraded Tevatron and in Fig. [5] for the LHC. (We note a disagreement 
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with Ref. ||12|| .) All calculations are again done at the leading order in QCD, and renor- 
malization and factorization scales are fixed as = fif = + m t . The cross sections 
are plagued by the same large uncertainties due to scale variations already observed in the 
case of the 2 — ► 3 processes. Notice that, away from the resonance region, m H ± < m t — mf,, 
the relative size of the two classes of cross section depends on yfs and m H ±. Indeed, at 
high energies, where the gluon-initiated 2 — > 3 processes dominate over the quark-initiated 
ones, the cross section arising from the elementary process gb — > tH~ is larger than that 
from 2^3 ones, which is penalized by a 3-body phase space suppression. At the Tevatron 
center-of-mass energy, the quark-initiated 2^3 processes still have the dominant role for 
values of m H ± that are not too large. For the particular choice of scales and fif made 
here, the cross-over for the two cross sections is at about m#± ~ 265 GeV. 



D. The inclusive single- production cross section 

Before presenting the inclusive cross section, some elaboration on the summation pro- 
cedure used to add the 2^3 and 2^2 channels is in order. Since the initial 6-quark is 
contained in the proton or antiproton via a gluon, the 2 — > 2 process is of the same order 
in a s as the 2 — > 3 ones. The collinearity of the 6-quark with the initial gluon induces 
the large factor a s (/x^) log(/^//mb), where the factorization scale [if is 0{m H ±) and it was 
chosen to be fif = m t + in our numerical evaluations. This factor is then resummed 
to all orders log(/i//mf>)) n when making use of the phenomenological 6-distribution 

function. The first order n — 1, is also contained in the set of 2 — > 3 partonic processes 
gg — > tH~b when one of the two initial gluons produces a pair bb collinear to the initial p or 
p. (See the last two diagrams of Fig. ^.) Thus, when summing the contributions from 2-^2 
and 2^3 cross sections, this term has to be properly subtracted in order to avoid double 
counting. Given the relevance of resummation for the large parameter a s (/i#) log(fif /rrif,), it 
is often concluded that the 2 — > 2 process gives the bulk of the cross-sections for the single 
production of H ± . As was already noticed in the previous sections, however, the issue of 
the dominance of one cross section over the other depends on yfs and m#± . 

The summation and subtraction procedure has been carefully systematized in the case of 



neutral Higgs bosons production |29| , |10| , for which, in a similar way, also the 2 — > 1 process, 
bb — > H°, and the 2^2 one, gb — > bH°, partially overlap. It was applied to the case of 
charged- Higgs production in Ref. ||30fl , where it was debated whether a 2 — > 1 elementary 



process tb — > H + also contributes to the inclusive cross section, once the theoretically 
calculated t-distribution function is folded in the proton beam. For the ranges of 
that may be probed at the Tevatron and at the LHC, the term log(m#/m t ) is sufficiently 
small, and the 2 — > 1 process can be safely omitted 0]. We therefore disregard the process 
tb — > H + in our analysis and in addition assume sufficient 6-tagging efficiency so that final 
states with three or more 6's can be distinguished P,[3l . 



Following Refs. pp| , pH| , we introduce the distribution function given by the per- 

turbative solution to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equation: 

i{x , flf)= ^Sl hi mf 1 ^P g J*-) g{ y, „,) (8) 

7T \m b J Jx y \y J 
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T q 




FIG. 7. The leading-order production cross section a(pp — * tH~X), for \fs = 2TeV, as a 
function of the charged Higgs mass, is shown for three different values of tan/3 = 2, 10, 50. The 
cross section is obtained by adding the contribution of the 2 — ► 2 process, gb — * tH~ , to those of 
the 2 — > 3 ones, gg — > tbH~ and — ► tbH~ , and subtracting overlapping terms. Renormalization 
and factorization scales are fixed as fiR = fif = m t + . 




250 500 750 1000 

m H - (GeV) 

FIG. 8. Same as in Fig. |^ for center-of-mass energy \fs = 14TeV. 
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where the splitting function is P g b(x) = (x 2 + (1 — x) 2 )/2 and g(x,fif) is the usual gluon- 
distribution function at the factorization scale /if. The hard process gb — > tH~ , convoluted 
with the distribution function b above, gives a contribution that has to be subtracted from 
the sum of the gg-initiated and (standard) gfr-initiated process, convoluted with the phe- 
nomenological 6-distribution function. 

The appropriately summed inclusive cross section is shown in Fig. |7| for the Tevatron 
and in Fig. |8| for the LHC As in the case of only the 2^3 channels, a potential reach 
of ~ 250 GeV and of 0(1) TeV is found, respectively, for the Tevatron (with 10-30 fb -1 
luminosity) and the LHC (with a luminosity of 100 fb" 1 per year). Background and other 
detection and identification considerations will somewhat diminish this reach ||, which is 
therefore given only as a rough guideline. 



III. SLEPTON-STRAHLUNG 



In supersymmetric models, lepton L and baryon B numbers are not accidental sym- 
metries, but have to be imposed by hand, e.g. by imposing a discrete i?-symmetry (im- 
parity) |32] Rp = (^_i) 3B + L + J ) where J is the particle spin. It is possible, however, that 
only L or B correspond to a conserved number, which is sufficient to ensure the proton 
stability. In particular, in supersymmetry the lepton and Higgs doublets are not distin- 
guished by their spin as in the SM. It is therefore natural to expect that some mixing exists 
between slepton and Higgs bosons as well as between leptons and higgsinos that carry the 
same quantum numbers, and, hence, that the lepton number L is generically not conserved. 
The realization of L-violation is basis-dependent and it is usually convenient to define the 
two-Higgs doublets as those whose neutral components are aligned along the two large vev's, 
whereas the lepton doublets are along the orthogonal directions in field space. In this basis 
it is straightforward to show that neutrino masses arise from small tree-level mixing with the 



neutralinos and at one- and two- loop levels from AL = 1, Yukawa-type interactions |33H36 
This offers an exciting avenue for generating neutrino mass and mixing at the weak scale. 
It further suggests collider tests of models of neutrino masses since both radiative neutrino 
masses and slepton production are controlled by the same Yukawa couplings. 

In the following, only models in which the lepton number is violated by AL = 1 renor- 
malizable operators are considered. The low-energy Lagrangian is derived from the super- 
potential operator 

W = X' lmn L l Q m D n , (9) 

where L, Q, D are the lepton and quark-doublet and down-singlet superfields, respectively. 
The possibility of renormalizable purely leptonic operators, also involving a AL = 1 lepton- 
number violation, does not affect our analysis, and such operators are neglected hereafter. 
In component fields, the slepton-quark interactions relevant to our purposes are: 

where / denotes the sfermion superpartner of a fermion / and all generation indices are 
arbitrary. In deriving (|10"D it was assumed that the right-handed quark fields as well as 
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the left-handed down-quark fields are already in the mass eigenstate basis, and that the 
CKM matrix V coincides with the rotation matrix of the left-handed up-quark sector. Had 
this assumption not been made, both fermions in each of the two terms of eq. ( JlOj ) would 
be multiplied by the appropriate rotation matrix. However, in the absence of any initial 
assumption on the texture of the A' matrix, this would merely correspond to a redefinition 
of its elements. 

The AL = 1 couplings given above cannot be arbitrarily large as they lead to tree-level 
corrections to various observables and correct neutral- and charged- current universality p7 



With the exception of A' m < 1CT 4 and A' lml ^ 10~ 2 from (/3/3) 0l ,-decay and atomic parity 
violation, respectively, one has X' Un < 0.01, X' l2n < 0.2, and X' l3n < 0.4-0.5, for sfermions 
with masses of 100 GeV. These constraints and their derivation are summarized, for instance, 
in Ref. []3~8|1 . As an example, the weak constraints on AJ 3n are derived from either t-quark 



decays [^,0 (this constraint, however, vanishes as the slepton mass approaches the t-quark 



mass), from the one standard deviation in the Z width [41], or from 6-quark semileptonic 
decays again at one standard deviation. The above constraints scale as a power of 

mj and are therefore significantly weaker for heavier sfermions. For example, for squarks 
near the 500 GeV mark, A' ~ 1 in the third family is generally not excluded. On the 
other hand, constraints on pairs of non-identical couplings (e.g. from meson mixing) often 
imply, further, that certain combinations of couplings cannot saturate their individual upper 
bounds simultaneously. In general, the hierarchy of couplings that emerges from experiment 
is similar to the generational hierarchy in the usual Yukawa couplings (with A 333 the most 
weakly constrained), an observation that we will adopt as a guideline. (A similar structure 
is also suggested by various theoretical models. See, for example, Refs. 



A. Production cross section 



The terms (|10"D in the Lagrangian lead to new and exciting possibilities for slepton produc- 



tion at hadron colliders: (i) exotic t-quark decays t — > rfe, if kinematically allowed |39| , |40| , |43 



(ii) s-channel resonant production of sleptons jH] (see, for instance, Ref. Q for a discussion 
of resonant production at LEP), and, as proposed above, (Hi) associate production qq'l. (In 
addition, (iv) gluonic couplings are induced for the sneutrino, which could now be singly 
produced gg — > v f45j.) Here, we focus on the charged-slepton-strahlung production, in par- 



ticular, stau t production in association with t- and 6-quarks. Associate production of the 



neutral sleptons will be discussed elsewhere [45]. Both 2 — > 2 and 2 — > 3 channels gb — > tf 



and gg, qq —> tf~b, will be considered (at the leading order) and properly combined. As in 
the case of the charged Higgs production, the 2 — > 3 processes encompass the production 
mechanism (i) in the relevant kinematic region. 

Concretely, the only .R-parity- violating coupling that is considered is X' 333 . This choice 
is motivated by the fact that this is the least constrained coupling, and, as explained above, 
it is expected to be the most significant among the AL = 1 couplings in some frameworks. 
In this case, the cross section scales as A333 in the kinematic region m t < nif + m^. In 
the complementary region m t ^ rrif + mj,, sleptons contribute to the width of the t-quark 
(in proportion to A 3 | 3 ), violating this simple scaling law. This implies that, although our 
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results may be taken as indicative for cases involving other A' couplings, they cannot be 
used directly in cases not involving the t-quark. In addition, the subtraction procedure to 
be followed when combining the different 2 — > 2 and 2 — > 3 channels into the inclusive 
cross section, although conceptually similar, differs technically for cases with one or two 
light quarks associated to the f-production. For this reason also, the results presented here 
cannot be simply adapted to production cross sections involving other A' couplings or to the 
case of neutral-slepton production. (The latter is also complicated by the presence of the 
gluon-fusion channel.) On the other hand, the slepton generation label does not enter our 
calculation but affects only the signal analysis (on which we comment below). Hence, our 
results can be generalized in a straightforward fashion to the production of other charged 
sleptons in association with t- and 6-quarks. 

It should be noted that the produced slepton is always left-handed, as a result of the 
structure of the operator (|H]). Of course, the physical eigenstates are, in general, admixtures 
of left- and right-handed current eigenstates. For simplicity, it is assumed in the following 
that the left-right mixing term in the f mass squared matrix is small. As a consequence, 
left- and right-handed current eigenstates coincide already with the two mass eigenstates, 
with masses vrif L and mf R . The first of these two masses, recurrent in this analysis, will be 
simply denoted as rrif. (A generalization is straightforward and involves the introduction of 
a mixing angle.) 

The inclusive production cross sections pp — ► tf^X, pp — > tfi,X are obtained by com- 
bining the production cross sections arising from the 2-^2 elementary process gb — > tf^ to 
those induced by 2 — ► 3 partonic processes, which give rise to pp — > tbf L X, pp — > tbf L X. The 
2^3 processes might be independently measured only if relatively complicated final states 
could be detected. The corresponding inclusive cross sections, formally given by eq. (fj), is 
obtained by convoluting the hard-scattering cross section of quark- and gluon-initiated pro- 
cesses with the quark and gluon distribution functions in p and p. The Feynman diagrams 
for the partonic processes are those of Figs. [I] and ^|, with H~ replaced by f^. Since the 
vector and axial coupling v and a for the vertex t-b-f^ are now simply v = a = 1/2, the 
square amplitudes |.M|jL and |-M|^ can be decomposed as: 

1 



M\ 2 q g=-\> 3 2 33 (Anas) 2 \V tb \ 2 (V« + A 



<re\ 



\M\ 2 gg = \\'£ 3 (Anas) 2 \V tb \ 2 (V« + A") . (11) 

The reduced square amplitudes V qq , A q ^ and V", A" coincide with those obtained for the 
2^3 charged-Higgs-production processes, once the replacement m#± — > m fL is made. 
In the kinematical region of a resonant t-quark, only the two decay modes t — > W + b and 
t — > f + b are considered. The latter has the width 

r(t-r + 6)=|2m, W (l + i-4)A^(l,4,4) , (12) 
327T \ vrq vrq J \ mi mi ) 

For simplicity, the charged Higgs H ± is assumed to be sufficiently heavy, as to kinematically 
forbid the decay mode t — > H + b. 
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If only the t-quark can be detected in addition to the decay products of ?l, the contribu- 
tion from the 2 — > 2 process has to be included. The Feynman diagrams for this process are 
those of Fig. || with the obvious substitution H — > f L ; the hard-scattering cross section 
is obtained from (|7D by replacing AGpM^/\/2 with A 333 /2 and m H - with m fL . The sub- 
traction procedure to avoid double counting follows exactly the pattern already described 
for the production of the charged Higgs boson. The final inclusive production cross sections 
a(pp,pp — > tfX) are shown in Figs. [5] and [TO] for the Tevatron and the LHC, respectively, 
as a function of the f L mass rrif L > 45 GeV. Note that, since lepton number is violated, we 
conservatively apply the model-independent lower limit rrif > 45 GeV extracted from the 
measurement of the Z width. 




100 200 
m~ (GeV) 



300 



FIG. 9. The leading-order production cross section a(pp — ► tfX), for y/s = 2 TeV, as a function 
of the f mass, is shown for different values of A 333 . Renormalization and factorization scales are 
fixed as iir = \i$ = m t + rrif- 



The large cross section obtained in the case of \/s = 14 TeV implies that at the LHC, 
with a luminosity of 100 fb _1 per year, light f's may be produced in abundance even for 
couplings as small as 0.01, whereas for large couplings, they may be produced up to masses 
of 0(1) TeV. 

The f-decay modes are highly model-dependent in the case of lepton-number violation. 
In particular, all superpartners typically decay in the collider and the typical large missing 
energy signature is replaced with multi-6 and lepton signatures, which may be used for iden- 

|l6f .) Depending on couplings and phase space, the main 

cb, ts, lu, where we included 



X°t, X ± 



u T , tb, 



tification. (See, for example, Ref. 
two-body decays for f are, for example: f 
the effect of purely leptonic couplings. (There also exist two-body decays due to tree-level 
Higgs-slepton, chargino-tau, and neutralino-neutrino mixing. However, these are strongly 
suppressed by the small mixing angle ~ {m v /mz) 2 '■) In addition, various three-body decay 
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F~t — i — i — r 



A 
I 



i — i — i — r=i 




250 500 750 
m~ (GeV) 



1000 



FIG. 10. The leading-order production cross section a(pp — > tfX), for y/s = 14TeV, as a 
iction of the f mass is shown for 
scales are fixed as fJ,R= Hj = m t + nij 



function of the f mass is shown for different values of A 333 . Renormalization and factorization 



channels may be open, depending on the model parameters. The charginos and neutralinos, 
if produced, also cascade in a mo del- dependent way to leptons and jets. Therefore, detec- 
tion and background studies cannot be done in a model- (coupling-) independent fashion; 
particularly so, once our assumption that only the charged Higgs or the stau (but not both) 
are produced in association with t is generalized. Nevertheless, many promising multilepton 
and multi-6 signatures are available. The study of detection aspects is well motivated by 
the potential reach in small coupling and/or large mass, but it is left for future works. 

Leaving detection issues aside, in L-violating models, sleptons are potentially more ac- 
cessible (depending on the coupling) than in models with lepton-number conservation, where 
their direct production relies on the Drell-Yan process, which allows a discovery reach for 
sleptons only up to mj < 350 GeV at the LHC [[47]]. (Cascade decays may provide the bulk of 
slepton production in the lepton-number-conserving models, but such processes are highly 
model-dependent. They also complement slepton production in the lepton-number- violating 
case studied here.) The potential reach at the Tevatron is limited to large couplings and/or 
light f's. Nevertheless, the slepton-strahlung provides a unique slepton-discovery mechanism 
at the Tevatron as well. 
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FIG. 11. A one- loop contribution to the Major ana neutrino mass arising from the AL = 1 
operators of the superpotential (||). 



B. Exploring AL = 1 models of neutrino masses 



An obvious question is to know to what extent does the discovery reach described above 
makes it possible to explore the corresponding models of neutrino masses. It is straight- 
forward to show that the supersymmetry-rotated operators ( |lOl) with a quark (sneutrino) 
replaced by a squark (neutrino) lead, at one-loop order, to a Majorana neutrino mass as 
illustrated in Fig. [□]. One obtains 

m v A , 2 /300GeV\ (v^\ ^ 



MeV y m~ b J ym&mj; 

where a 5-quark and 6-squark are assumed to circulate in the loop, and m 2 LR is the b left-right 
mixing squared mass. (Note that the size m 2 LR here can be significant even if left-right stau 
mixing is suppressed.) The neutrino mass may vanish in the limit of a continuous U(l)# 
symmetry, which corresponds in our case to m\ R <C m^m^. It is further assumed that 
no other accidental cancellation among various contributions to the neutrino mass takes 
place. In this case, the contribution ([13]), if it exists, constitutes probably the dominant 
contribution to the neutrino mass (another contribution may arise from tree-level neutrino- 
neutralino mixing). Imposing laboratory limits on the v e mass one can, for example, derive 
severe constraints on A' 133 [ §j| . 

Given all other constraints and the cross sections of Figs. |9| and |1(], A233 and A 333 are 
the couplings that are likely to be probed at hadron colliders through slepton-strahlung 
production. Assuming that associate stau production can be observed for A' ^ 0.1 (0.01), 
neutrino masses heavier than 3 KeV (30 eV) can be explored, for < 1 TeV and mj mg/ m 2 LR 
of 0(1) in (|I3"D . It is interesting to note that if the coupling is large enough to lead to slepton 
production, it is probable that the corresponding neutrino cannot be in the sub-eV range, as 
is sometimes assumed, unless the 6-squarks are in the multi-TeV range, and with negligible 
left-right mixing. Furthermore, a large slepton production cross section, may even imply 
a heavy neutrino species, which decays on cosmological time scales. Alternatively, if the b- 
squarks are not discovered at the LHC (when considering the -R-parity- violating cascades), 
the same range of coupling would automatically imply much lighter neutrinos, possibly in 
the sub-eV range, enhancing the coverage of the neutrino-mass parameter space. 

Thus, collider studies in this case carry indirect but crucial implications to models of 
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neutrino masses and can help reveal the neutrino spectrum. This beneficial relation is due to 
super symmetry, which relates the neutrino radiative mass operators and the slepton-quark 
Yukawa operators. Negative search results can alternatively provide strong constraints on 
the A^ 33 couplings, especially if slepton masses are independently measured. Such potentially 
strong constraints are currently not available by any other method. Further handles on the 
couplings and on the neutrino spectrum are provided by the search for the sneutrinos [ 45j , in 



which case also the gluon fusion gg — > v channel is available. The gluon-fusion cross section 
also depends quadratically on A' Yukawa couplings, and it is given by a straightforward 
generalization of the Higgs gluon-fusion gg — > H° cross section [19]. 



IV. SUMMARY AND OUTLOOK 



In summary, we have shown that charged Higgs- and slepton-strahlung provide important 
channels of single production. Our study implies a significant production cross section, in 
particular at the LHC, but also at the upgraded Tevatron energies. For example, charged 
Higgs bosons and staus as heavy as 1 TeV can be produced at the LHC (Figs. ^ and |1(]). For 
lighter staus, .R-parity- violating couplings as small as 0.01 may be probed (Fig. |H]). At the 
Tevatron (Figs. [7] and [|), the kinematic reach may be significantly extended with respect to 
that obtained from the t-quark decays t — > H + b and t — > f + b. Of course, more conclusive 
statements should await detailed background and detection studies. 

Our calculations, at the leading order, include exact evaluation of the three-body phase 
space, and special attention was given to the correct summation of the various contributions 
to the inclusive cross section. 

The importance of 6-tagging in separating the various channels was pointed out. In 
addition, various weakly interacting particles may be produced simultaneously in association 
with quarks (e.g. the charged Higgs and one or two charged sleptons), hence raising the issue 
of their separation and identification. In particular, in the case of the stau studied here, 
the slepton signature could be similar to that of the charged Higgs. (Similarly, the neutral 
Higgs and sneutrino could both be produced via gluon fusion and/or strahlung and could 
also decay similarly, extending the problem to the neutral sector as well.) One may have to 
rely on mass differences and on more suppressed L-violating or supersymmetric decays in 
order to distinguish the different bosons. 

Detection of singly produced sleptons via the strahlung process (or any other process) 
carries substantial benefits to the mapping of the lepton-number-violating potential (and 
superpotential). Hence, it also carries important consequences to models in which neutrino 
masses are obtained radiatively, since radiative neutrino masses are proportional to the 
couplings of the slepton-quark Yukawa operators. Leaving detection issues aside, heavy 
sleptons could be abundantly produced in L- violating models and small L- violating couplings 
(and hence, small neutrino masses) may be probed. 

Though not discussed explicitly, a similar situation to the one studied in this paper can 
arise in any other model in which a weakly interacting scalar couples via a non-negligible 
Yukawa coupling to quarks. The most obvious example, which corresponds to a straightfor- 
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ward generalization of the slepton-strahlung case, is given by (scalar) lepto-quark models. 
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APPENDIX A: THREE-BODY PHASE SPACE 



In eq. (Q) dPS (q± + q 2 ; Pi, p 2 , P3) is an element of the 3-body final-state phase space 
normalized as ITi 



dPS 



(2vr) 5 (qi + q 2 -Pi-P2-P3) 



1 1 r +1 



d 3 pi d 3 p 2 d 3 p 3 



(2vr) 3 2Ei (2vr) 3 2£ 2 (2tt) 3 2E 3 



r2ir fi?i,max i-i?i2,max 

, d(cos9) / d(f) dEi \ dE 12 

(2tt ) Si-i Jo Je l,min J_Ei2,min 



(Al) 



where qi (p^ are the four- momenta of the initial partons (final-state particles) and Ei, E 2 
can be chosen, e.g. as the energies of the final t and b respectively, with E\ 2 = (Ei — E 2 )/2, 
and E3 the energy of the H . The remaining integral in eq. (SI) is over appropriately 
defined angles 9 and 0, describing the motion with respect to the beam axis of the three- 
momenta pi, P2 of the two produced quarks. (Of the initial four angular integrations, one is 
eliminated from energy conservation and another (azimuthal) angle integration simply gives 
a factor of 2ir included in eq. (|A1|) .) The integration bounds for the energy E\ are as follows: 



s + ml — (m b + m H ±Y 
2V§ 



(A2) 



those for the energy E 2 : 



E 



12,min/max 



-b =F Vb 2 - 4:ac 
2a 



(A3) 



with a, b, and c given by: 

a = 2EiVJ — s — m\ 



\fl - Ei) (m 2 H ± - m\ 



m. 



VI -E x ) 2 -Am 2 b ~ (El + m 2 H± 



m. 



m, 



(A4) 



Finally, the resulting six-dimensional integral eq. (|J) over the remaining phase space and 
over the parton luminosities is performed numerically with the standard Vegas Monte-Carlo 



integration routine 49 . 
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